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Abstract 



o i 

A lattice action for QED is considered, where the derivatives in the Dirac 



operator are replaced by one-sided lattice differences. A systematic expan- 
sion in the lattice spacing of the one-loop contribution to the fermion self 

energy, vacuum polarization tensor and vertex function is carried out for an 
Q-i' 

arbitrary choice of one-sided lattice differences. It is shown that only the 



vacuum polarization tensor possesses the correct continuum limit, while 



^ . the fermion self energy and vertex function receive non-covariant contri- 

butions. A lattice action, discretized with a fixed choice of one-sided 
lattice differences therefore does not define a renormalizable field theory. 
The non-covariant contributions can however be eliminated by averaging 
the expressions over all possible choices of one-sided lattice differences. 
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1 Introduction 



As is well known, the naive discretization of the fermionic action which preserves 
the antihermiticity of the derivatives in the continuum, leads to the notorious fermion 
doubling problem. One of the most popular proposals for eliminating the fermion- 
doubling problem consists in working with Wilson fermions Q, whereby the action is 
modified by a so-called irrelevant term which breaks explicitely the chiral symmetry. 
This is consistent with the Nielson-Ninomiya theorem which states that for an action 
respecting the usual translational invariance, locality and hermiticity requirements, 
the fermion doubling problem can only be eliminated by breaking chiral symmetry . 
Another way to circumvent the above theorem is to give up reflection positivity. Thus 
replacing the derivatives in the continuum by one-sided lattice differences ensures 
that high momentum excitations arising from the corners of the Brillouine zone do 
not contribute in the continuum limit. In fact the functional formalism of Berezin || 
dictates that the temporal derivative in the Dirac operator should be discretized using 
the one-sided lattice differences. By symmetry one would then be led to a similar 
discretization of the spatial derivatives. In this case the Hamiltonian would no longer 
be hermitian for finite lattice spacing, which does, however, not necessary imply that 
correlation functions do not possess the correct continuum limit. Actions involving 
one-sided lattice differences have been considered in [Q, |5], |J, where by a suitable 
averaging procedure, reflection symmetry is effectively restored a posteriori. 

In this paper we will study in detail the behaviour of the fermion self energy, 
vacuum polarization tensor and vertex function for small lattice spacing in one-loop 
order, for a lattice action discretized with one-sided differences. This analysis is car- 
ried out using the small-a-expansion scheme discussed in ]7). Such a systematic study 
has not been presented in the literature, which is surprising in view of the importance 
of the problem. It is shown that while the vacuum polarization has the expected 
continuum form for small lattice spacing, the fermion self energy and vertex function 
include non-covariant contributions which can only be eliminated by averaging the 
expressions over all possible choices of one-sided lattice differences. These compu- 



2 



tations are carried out in sections |^ and Since the small-a-expansion scheme of 
ref. |7j has not been widely used in the literature, we shall discuss it in detail in the 
following section. This section also serves to define our notations used in sections |] 
and |. 

2 Small— a— expansion 

In this section, we discuss in detail the small-a-expansion scheme of ref. |7[] 
following closely the proof given in that reference. 

A lattice regulated one-loop Feynman integral in D dimension has the form: 

+^ 

f d D p' 

a(p,m,a,D) = J D H (p',p,m,a,D) , (2.1) 

a 

where p stands for the set of external momenta, p' is the loop momentum and a the 
lattice spacing. The particle masses are denoted collectively by m. An important 
feature of the lattice regulated Feynman integrals is that the lattice spacing appears 
in the integrand as well as in the integration limits. This integrand is a periodic 
function of p and p' and possesses a continuum limit. Furthermore H (p',p, m, a, D) 
has the following homogeinity property, 

H (p',p,m,a, D) = a dn H (ap' , ap, am, a = 1, D) , (2.2) 

where du is the inverse-mass dimension of the integrand. In the naive continuum 
limit the integral (|2.1|) reduces to the usual, in general divergent, continuum Feynman 
integral. Since H (p',p,m,a,D) is an analytic function of a in the neighborhood of 
a = 0, we can expand it into a Taylor series. Consider the expansion up to O (a J ) , 

j 

H (p',p, m, a, D) = ^2 aJ Hj (p',P, m ; D) + Rj (p',p, m, a, D) . (2.3) 

i=o 

The Taylor coefficients Hj (p' ,p,m, D), which have a structure resembling that of 
the integrands of continuum Feynman integrals, possess the following homogeinity 
property: 

Hj(ai/,ap,am,D) = a j ~ dH Hj {p' ,p,m, D) . (2.4) 
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From ( |2.2| ) and ( |2.4| ) it further follows for the remainder Rj that: 

Rj (p ,p,m,a, D) = a d " Rj (ap ,ap,am,a = 1, D) . (2.5) 

Introducing ( |2.3|) into we have that: 

J + f d D p' 

a(p,m,a,D) = ^ a 5 (p, m, a, D) + / D Rj {p',p,m,a,D) , (2.6a 

j=0 A ( 27F ) 

where 



f dPp' 

(?j{p,m,a,D) = a J J D H 3 (p',p,m,D) . (2.6b) 



a 

Consider first (|2.6b|) . This expression can be decomposed as follows: 



<jj (p, m, a, D) 



+oo 

-.0 



f d D p' . f d D p' . 

-oo ^ ' |p'|>- 



(2.7) 

For j = 0, the first term corresponds to the naive continuum limit of ( |2.1| ). For D = 4 
this integral in general diverges. Hence from now on it will be understood that all 
integrals are regulated using the dimensional regularization schemeQ Notice that for 
a / and m ^ the original lattice Feynman integral ( [2.1| ) will be well defined. 
Hence any divergencies introduced by this decomposition have to cancel eventually. 

Making use of the homogeneity property ( |2.4j) , and introducing the dimensionless 
integration variable p' = ap f , expression ( |2.7| ) takes the form: 

<jj (p, m, a, D) = 

dP p ! f dF* p ! 

= a J ~F?~F Hj p ' m ' D ^ ~ adH D J Hj aPl am ' D ^ ' 

-oo ^ 7 1 > |p'| >7r V n ) 

We next expand i/j (p', ap, am, D) in the lattice spacing a up to O (a J ~^ dH ~ A ^. Then 

<7j (p, m, a, D) = 

/~^~ dF* p ! f d^ p f 

(a J+ ( 4 - D ) +1 ) , (2.8) 
*Dimensional regularization of lattice Feynman integrals have also been considered 
by Kawai et al. M 



where Tkf (a) stands for the Taylor expansion of / (a) around a = up to order k. 
The expression ( |2.8|) can be trivially rewritten as follows 



<7j (p, m, a, D) 



' d D p' 



Hj{p',p,m,D) 



a 



d H -D 



d D p' 



Tj-(d H -A)Hj (p', ap, am, D) 



-a 



d H -D 



J j^s Tj-te-QHj (jf, ap, am, D) + O (a J+ ^ +1 ) . (2.9) 



The p'-dependence of the coefficients in the Taylor expansion of Tj_^ H -4)Hj have the 
form p'^ • • -p' H l (p' 2 ) ■ Hence the last integral vanishes in the dimensional regulari- 
zation scheme. Introducing ( |2.9| ) into (|2.6a|) , we therefore have that: 



a (p, m, a, D) 



+00 



a-' 



+a 



j 

E 

i=o _ x 
+- 

d H -D 



d D p' „ _ _ ^ , _ dB - D 7 d °P' 



Hj {p',p,m,D) + a a 



q Rj (p', ap, am, a = 1, D) + O (a J+{4 " D)+1 



■ d D p' 



(27T) 



D Tj_ (dH _ 4) if j (p', ap, am, L> 

(2.10) 



where we have made use of the homogeneity property ( |2.5| ), and have introduced the 
dimensionless loop momentum p' = ap' in the last integral in ( |2.6a ). From 
and (12.51) it follows that: 



Hj (p, ap, am, D) = H (p', ap, am, a — 1, D) — Rj (p', ap, am, a — 1, D) 

3=0 



Introducing this expression into ( [2.10 ), we therefore have that: 



J ■ 7° d D p' 
a (p, m, a,D) = ^a 3 / — Hj (p', p, m, D) 
3=0 (2vr) 



+a 



4-D 



a 



Hence the small-a-expansion scheme can be summarized as follows: 



a (p, m,a,D = 4) 



lim 



(p, m, a, D) + a 4 u a (p, m, a, D 



4-D~ 



(2.11a) 
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where 

j 



with 



<7 (oo > (p, m, a, D) = ^aV (oo)/(i) (p,m, J D) + o(a J+1 ) 

3=0 



„M/0, (p , m , g) = 1_JJL^ ( ggg^i ) , (2 .„ b ) 



a=0 



and 



where 

(p, m, D) 



j 

a{p,m,a,D) = £ a j d (j) (p, m, D) + C (a J+1 ) 

j=d H -4 



1 + / d^p' (d j - dH+4 H (p', ap, am, a = 1, £>)* 



(j-d H + A)\J (2tt) d V da^+ 4 



I a=0 

(2.11c) 



3 Small— a— expansion of the fermion self energy 



The action of the lattice U (1) gauge theory we shall consider is taken to be of 
the form 

/' 

(/ 



Sty, if); A; {e,}} = S G [A] + £ if) (x) \m--J2 Vfr - 1 



11=1 



if){x), (3.1) 



where I] = I] a 4 , and we have replaced the derivatives <9 M in the fermionic contribution 

x n 

to the continuum action by one-sided lattice differences || according to: 

d^ip(x) -> -I if) (x + ae^p) - if> (x) j = -e M (T eiitJ , - l\ if) (x) . 

Here /a is a unit vector pointing in the /^-direction, and e M = 1 (e M = — 1) corresponds 
to choosing the right (left) lattice difference. For = — 1 the link-variable LL^ is 
defined by (a?) = UAx — ap). The contribution Sg is the standard plaquette 
action for the gauge field. From (|3.1| ) one readily obtains the following expression for 
the fermion propagator, the photon propagator in the Feynman gauge, and vertices, 
coupling one or more gauge potentials to a ip-ip pair: 

6 



Fermion-Propagator: 



S F (a,p,m, {e M }) = 
Photon-Propagator: 

Vertaces: 



4 

° M=l 

4 

m 2 + 4E sin 2 ^e ie ^" a 

p=l 



(3.2a) 



(3.2b) 



V^(a,p, {e M }) 
(a,p, {e M }) 



+ia 2 g 3 5 fl2lll 5 lllfllfl e ie ^ a . 



(3.2c) 



Here P M = 2 (P m + P° ut )^ where p^ 1 , is the incoming (outgoing) momentum of 

the fermion. 

The fermion self energy receives a contribution from the two diagrams depicted 
in fig. |l|. The corresponding dimensionally regulated Feynman integrals are given by 



d D p' 



S (ct) (p,m,a, {e^},D) = #V D / ~~~d (p' , p, m, a, {e^}, D) , (3.3a) 



where a = a,b and 



ffW (p',p,m,a,{e M },L>) 



E 



4 £ • 

? E sin 

T=l 



2 (p-p') T a 
2 



m 



2 + Jr E sin 2 ^ e ^ a 



(3.3b) 



and 



D 



H^(p',p,m,a,{e,},D) = -a £ e^e lw 



£ E sin^ ^ 



(3.3c) 
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In ( |3.3a| ), /i is the usual arbitrary mass scale introduced in the dimensional regulari- 
zation scheme. 

According to fl2.11a| ) we have: 



£(<x) (p,m,a, {e^},D = 4) 



lim 



(3.4a) 



where and S^j are defined by 



<?y- D / ^yn OA P, ™, o, M, D) + (a J+1 ) , (3.4b) 



and 



£ ((j) (p,m,a, {e M },£>) 

1 2 7 



^ 2 / ^+i^ (,T) (p', aP, am, a = 1, {e M }, £>) + (a J+1 ) . (3.4c) 



a 

— 7T 



(27T) 

Here T[H^ denotes the Taylor expansion of H a in the lattice spacing a around a = 
up to order O (af^j . All integrals are to be calculated using dimensional regulariza- 
tion. In the limit D — > 4 the coefficients of a J in the expansion of ( |3.4b|) and ( |3.4c|) 
will, in general, be ultraviolet (UV) and infrared (IR) divergent, respectively. The 
UV divergencies can be isolated in the standard way, since the integrals are of the 
continuum type. The technique for isolating the IR divergencies is described in the 
appendix. For o ^ 0, the UV and IR divergencies must cancel, since the original 
lattice Feynman integrals are finite. 

Consider first the leading contribution for a — > to (|3.4a ), which is deter- 



mined by the coefficient of O (a x ) of ( |3.4c| ). This coefficient, which we denote by 



({e M },D), receives contributions from both diagrams in fig. |l|. After making the 
change of integration variables e^p'^ — > p' , one is led to the following non-covariant 
expression 

D 

({e M }, D) = J2 W % (D) , (3.5a) 
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where 



D 



a, (D) 



-9 



7 d Dp, 



4 E sin z y -f 
P =i 



D 

E 



+ 1 



4 E sin 2 v -fe ip 'p 



(3.5b) 



is IR-convergent. Note that by averaging (|3.5a|) over all possible 2 D sets of {e M } this 
expression is seen to vanish, as has already been pointed in ref. ||. 

Next consider the contribution of O (a ). In this order reduces to the usual 
continuum form of the Feynman integral for diagram (a) in fig. [I], while Ej; 00 '* does 
not contribute in this order, as follows from ( |3.3c| ). Hence £(°°)/( ) is given by the well 
known dimensionally regulated continuum expression for the fermion self energy, 



S (°°)/(o)( p ) = A {co) (p) + (m - ift (p) , 



(3.6a) 



where 



A {co) (p) 



B {oo) (p) 



+ 



3mg 2 



?tt 2 (4 - D) 16tt 2 



, m 

H 2 



mg- 



± 

J da (1 + of) In 



am 2 + a (1 — a) p 2 



a 2 m 2 



+ 



q 2 , m 2 
In — 



8vr 2 (4-D) 16tt 2 /i 2 



da (1 — a) In 



am 2 + a (1 — a) p 2 



a 2 m 2 



-|- • 



(3.6b) 



(3.6c) 



The constants C^°°^ and Cg°^ are given by 



C 



(bo) _ 3m# 5 



167T 2 

with 7^ the Euler constant. 



7s - 



5mg 2 

167T 2 



✓7(00) 



16n 



-,1E 



16tt 2 



(3.6d) 



Consider next the contribution of O (a ) to £ = S a + £&, where E( CT ) (cr = a, 6) 
has been defined in ( |3.4c|) with H^ a ' given in (|3.3b| , c). Both diagrams (a) and (b) 
contribute in this order. We denote this contribution by £(°). It is given by 



t^(p,m,{e,},D) = ^ ) (p,m,{e,},D) + tf\p,m,{e,},D), (3.7a) 



(0) 



where 



Si' 



(o) 



'2g 



f aD ~, E e^sin^e P 2 M [2e M p /1 e < ^i' - £ e v p v e lp> 

2 [ C( J) /i=l \ ^=1 



(27T) 



1? 



4 E sin 2 ^ 
. P =i 

1 



* d D p' 



D 



4 E sin 2 ^ e ip 'p 
. P =i 



(2 



7T 



E 

P =i 



4 E sin 2 ^ 



(3.7b) 



is infrared convergent for D = 4, and 



+ 7T 



-mg 



d D p' 



d .., 



(2 



7T 



D 



4 E sin 2 



P =i 



D 



4 E sin 2 v -f e ip '. 



D 



p=l 
D 



D 



-Aig 



2 / d D p' A=l 



, E e\P\sinp\ E e^sin-^e 2 2e jp f - E e 

p,=l \ u=l 



IP , 



(27T 



D 



4 E sin 2 ^ 

. P =i 



4 E sin 2 v -f e ip 'p 



(3.7c) 



is infrared divergent for D — > 4. The second term on the rhs of ( |3.7b| ) is the contri- 
bution from the diagram (b) in fig. [I|. After some manipulations, expressions ( ft.7b|) 
and (|3.7c ) can be written in the form: 



Ef (p,m,{e M },L>) = A t + (m - if) B t + C % (p, {e M }, D) 



1,2. (3.8) 



Here Ai and i5i are e^-independent finite constants, involving lattice integrals which 
can only be computed numerically. The e^-independent coefficients A 2 and £> 2 are 
infrared divergent for D — > 4 and have the form 

A 2 = -3mg 2 M (D) + V A , 
B 2 = - g *M(D)+V B , 



where T>a and V B are finite constants, and M. (D) is given by (see eq. ( |A.5| ) in the 
appendix), 

d D p' 



M(D) 



1 



1 



1 



(27T 



D 

E 



4 E sin 2 ^ 



8tt 2 (4-D) 



■ A2 



(3.9) 
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with Mi a numerical constant. 

The e M -dependent part in (|3~8| ) has the following non-covariant structure: 
Ci(p,{e„},£>) = 

d D P ' 



p+V 

C 2 (p,{e^},D) 



pi l V n ■ 

sin e 2 e jp w 



(2vr 



D 

4 X) sin 

P =i 



2 P_p 
2 



D 



4 X sin 2 ^ e ip 'p 
P =i 



(3.10a) 



pi ip' u I ■', D ,- / 

jr D ~ f smp / u sin h^-e 2 2e tp f — £ e ip 



dp 



A=l 



(27T) 



4 £ sin 2 ^ 



4 £ sin 2 ^ e ip, < 



L p=l 



-.(3.10b) 



While ( |3.10a| ) is finite for D = 4, ( p.lUbp is infrared divergent for D — > 4. 

From (|3.7a|) - (|3.10b|) we therefore find that, in this order the contribution of the second 

term on the rhs of ( |3.4a| ) to £ = S( a ) + £(&) is given by: 



i+(m-^)B + C(p,{e /i },D), (3.11a) 



where 



.4 
B 



3mg 2 1 3mg 2 . ~ 
h ^r^- In {ail) + C A , 



8tt 2 (4 - D 



■in' 



In (a/i) + Cb- 



(3.11b) 
(3.11c) 



2 _, 



8tt 2 (4-D) 8vr 2 

Here Ca and C# are finite constants, and C (p, {e^}, -D) = £ Cj (p, {e^}, -D) (see eqs. 
i=i 

(|3~T0a| , b)). Note that the expressions (|3.11b| , c) involve two regulators: the lattice 
spacing a and the dimension D. 

Combining the contributions p.5a| ), ( |3.6a| ) and ( |3.11a| ) we arrive the following expres- 
sion for the fermion self energy for D = 4, 



S(p,m,a, {e M },£> = 4) = 
^ 4 4 
= - H + tp.tvlnPifrp.v + A+{m-ip')B + 0(a) 

a H=l n,u=l 

p^u 



(3.12a) 



where a M = a M (Z) = 4) and the (for D — > 4 infrared divergent) coefficients are 
obtained by combining ( |3.10a| ) and (|3.10b| ). The expressions for ^4 and B are given 

by, 
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A 



B 



8m + 



rag 



sir' 



da (1 + a) In 



am 2 + a (1 — a) p 2 
a 2 m 2 



1 - 



SIT' 



da (1 — a) In 



am + a(l — a)p 



a 2 m 2 



(3.12b) 



where 



8m 
ZT 1 



3mg 2 



In (am) + 



4 Z 



■ In (am) + C 2 . 



(3.12c) 



Here C m = Cjf + C A and C z = eg* + C B , are the e^-independent constants. These 
constants differ from those obtained in the dimensional regularization scheme for 
continuum Feynman integrals, since we have used the lattice clS db TG gulator. 

Note that a remarkable cancellation of e^-independent, but in the limit D to 
4, infrared and ultraviolet divergent terms has occured after combining the above 
mentioned contribution, leaving us with a /i-independent expression with the lattice 
spacing as the only regulator! 

The first two terms in ( |3.12a| ) cannot be eliminated by introducing a renor- 
malized mass parameter, and renormalized fields. Both terms have a non-covariant 
structure for any choice of {e^}. These non-covariant contributions are seen to va- 
nish if they are averaged over all possible sets of {e^}. Such an averaging procedure 
had been proposed in ref. ||, where the fermion self energy has been discussed on a 
qualitative level. 



4 Small— a— expansion of the Vacuum polarization 
tensor and vertex function 



We begin with the discussion of the small-a-expansion of the one loop contri- 
bution to the vacuum polarization tensor. 
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i ) Vacuum polarization tensor. 

The diagrams contributing to the vacuum polarization are shown in fig. |2], of 
which diagram (b) has no continuum analog. The corresponding Feynman integrals 
have the form: 

+- 

U^(p,m,a,{e,},D) = g 2 ^~ D f —^H^(p',p,m,a,{e^D), a = a,b, 

_2T 
a 

(4.1) 



where the integrand w can be readily obtained using the expressions for the 
propagator and vertices given in ( |3.2a| -c). The small-a-expansion of IT^ analogous 
to ( |3.4a| -c) now reads: 



n^ ) (p,m,a J {eJ,D = 4) 



lim 



n (a) (P, m , a , { e ^}> D) + (a/i) 4 D IL^ (p, m, a, {e^}, D 



(4.2a) 



where and nj^ are defined by 

lC/ (oo) (p ) m,a ) {c M },£>) = 



+ 00 j£) , 

9 2 ^'° J ^ 5 T J H^(p',p,m,a,{e ll },D) + 0(a J+1 ), (4.2b) 



and 



n fa) (P, ™, a, {e^}, D) 



a 



\g 2 J -^n Tj +2 H»»/W (p', ap, am, a = 1, {e,}, D) + O (a J+1 ) . (4.2c) 



> y2 J (2VT 

— 7T V 



The leading contributions for a — > to ( 4.2a| ) is determined by the coefficient of 



O (a 2 ) and (9 (a x ) of which we denote by 2 ^ and ft^/ ^ ■ The coefficient 



of fi/^ 2 ^ for o = a and a = 6 is given by 



nf /( " 2) ({^}^) 



_7r ,n", • p'„ • b' S • */ 3 ■ *i 

f d p' sm sm sfe2 ip ^ea ip " 
-32g%e v ' - 



( 27r ) D [4Esin 2 ^e^ 



(27f 



P =i 



4 £ sin 2 ^e ip 't 
p=i 



(4.3a) 
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and 



n 



W(~2) 



D) 



Ug 2 5, 



d D p' 



■ P'u 

sin -£-ei %p * 



D 



4 £ sin 2 ^e ip 'f 



(4.3b) 



After a partial integration of the first term in ([4.3a ) one finds that 



fW(-2) = _fx 



W(-2) 



(4.4) 



Hence there is no quadratic divergence (for a — > 0), as expected from gauge invariance. 
The lattice expression for the coefficients ft^/ is rather lengthy and we do not 
present it here. After some trigonometric manipulations and partial integrations we 
find that 



■nW(-l) 



ttM^/(-1) 
U (6) 



0. 



(4.5) 



Hence there is also no linearly divergent contribution to U' 1U for a — > 0. Note that 
(|4.4| ) and (|4.5|) hold for arbitrary choice of {e^}. 

Consider next the contribution of O (a ) to H^ / /^°°) in eq. ( 4.2a ). For a = a it 
is given by the usual dimensionally regulated continuum expression for the vacuum 
polarization tensor, while n{^°°^ does not contribute: 



r /W ]o(rf) = -PV) n(oo)/(0) (P,™,D) 



(4.6a) 



where 



n (oo)/(o) ( P)m)jD ) 



1 



6tt 2 (4-D) 12tt 2 



In 



m 



2vr 2 



a (1 — a) In 



m 2 + a (1 — a)p 2 



da + C {QO \ 



(4.6b) 



and C(°°) 



12^ 7£- 

The expression for the contribution of O (a ) to H.^ 



n 



l (a)^ n (fc) is very lengthy. 
After some work (involving mainly partial integrations) one finds 

Pa / 7^ 

A=l 



n^/(o) ( P) D j 



(2^ 



4 £ sin 2 ^e ip 'p 

3 2ip / /lp 2ip / „ 



4 £ sin 2 ^e ip 't 
. p=i 



(4.7) 
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This contribution can be shown to be of the form 



nw(o) (Pj D) = ^ Pu _ p 2 S ^ n(o) (jD) 



(4.8a) 



To this effect one must consider separately the cases when the momenta appearing in 
the exponentials of the integrands carry the same or different indices. Making use of 
the fact that for p ^ A the integrals in ( |4.7|) are independent of the choice of p and 
A, one finds that 



m (d) 



-g 2 M (D) + 



fiW (D) 



reg. 



(4.8b) 



where A4 (D) is given by ( |3.9|) . The second term in ( |4.8b| ) is regular in the limit 
D — > 4, and reads 



ri (0) (d = 4) 



reg. 

d 4 p' 
(^j 1 



12 



E 7 



1 



/i , (7 = 1 



4 £ sin 2 ^e ip 'p 
L P =i 



4 £ sin 2 ^ 

. p=i 



(4.8c) 



We therefore find that the contribution of the second term in 



in this order, is 



given by 

M 4 - D rp^°) (p, d) D « 4 ( PmP „ - p 2 5 



-s 2 i s 2 



In (ap) + C 



(4.9) 



,6tt 2 (4-D) 6tt 2 

where C is a finite constant. Notice again that this expression involves two regulators: 
the lattice spacing and the dimension D. By combining (|4.9| ) with ( |4.6a| ) the in- 
dependent terms are seen to cancel, and we are left with a p-independent expression 
valid up to order O (a ), in which the lattice spacing appears as the only regularization 
parameter: 



n"" = (p»Pv - P %„) n<°> (p) + o (a) . 



(4.10a) 



Here 



n<°> (p) 



Z3- 1 - 1 



2tt 2 



a (1 — a) In 



m 2 + a (1 — a)p 2 



da, (4.10b) 
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and 

a 2 

Z 3 X -1 = In (am) +C, (4.10c) 

where C = + C is a finite e^-independent constant. Notice that in contrast 
to the case of the fermion self energy (|3.12a|), where the non-covariant contributions 



of O (a -1 ) and O (a ) could only be eliminated by averaging the expression over all 
possible sets of {e^}, the expressions ( |4.10a| -c) have the correct continuum structure 
for any choice of {e^}. The expected cancellation of the D-dependent terms, observed 
above, has been checked up to O (a 2 )f\. 

ii ) Vertex function. 

In analogy to ( |3.2c|) we write the vertex function in the form: 



«E/J(P + P') M 

V^ 1 ' {p,p',m;a, {e M },£>) = —ig (7^ + A M (p,p r , m; a, {e /t }, D)) e 2 , 

In one-loop order A M receives a contribution from the diagrams shown in fig. [3|. The 
lattice Feynman integrals contributing to A M , which we denote by (er = a, b, c, d) 
can be readily written down using the expression ( |3.2a| -c) for the propagators and 
vertices. The small-a-expansion of A M = E now reads: 



A M (p,p',m;a, {e^.D = 4) 



lim 

£>-*4 



(4.11) 



with Ajf ) = EA|T )/W and A M = EAj^. Here Ajf5/(o°) an d Aj^ are defined by 
expressions analogous to ( |3.4b| ) and ( |3.4cj ) , respectively, except that the factor a -1 



in ( 3.4c|) is replaced by 1, and Tj +1 is replaced by Tj. The leading term in their 
expansions in a is of order O (a ). In this order A™'°°' is given by the usual D- 
dimensionally regulated continuum expression, while K^'^°°^ vanishes for a = b,c,d 



in this order. Hence 



f Only one of us (N.S.) had the nerve to carry out this very extensive computation. 
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where A^ e ^ is the known regular contribution of the continuum formulation. 

Next consider the contribution of O (a ) to A M . All four diagrams in fig. |] 
contribute in this case. The expression for Aj^ is rather lengthy, and we will not 
give it here explicitely, since it can be obtained in a straightforward way. After some 
simple algebra, one find that A^ has the following momentum independent form: 

A? 1 = ^7, + 7/.£ + C ({£,},£>). (4-13) 

The first term is entirely determined by diagram (a). The second term, which receives 
contributions from all four diagrams, is finite for D — > 4. Finally, the third term, 
appearing on the rhs of (|4.13|) , is a non-covariant, e^-dependent term. It is given by 



D 

+ A 9 2 J2 e 8 e rlvlrl^lelv J 



(jP^, e'""e""' l sm^sm=e 2 1 



D 



(2vr 







YD t 


2 


4 £ sin 2 ^ 




4 £ sin 2 ^-e ik p 




L P =x 




I P =i _ 





d D k 
7d 



sin^e 2 e ^ 



r=l 









4 £ sin 2 ^ 




4 £ sin 2 %e**" 


L P =i 




L p=i J 



(4.14) 



where the first term arises from diagram (a), and the second term from (b) and (c). 
This expression can be reduced to the form 



D 



D 



Cp ({e M }, D) = J2 ^rlfTlrle (inure (D) + J2 e ^ e rlr V (D) , (4.15) 



V,T,6=1 

r^G 



T=l 



where the coefficients a^ UT g (D) are given, in the limit D — > 4, by infrared divergent 
integrals, while b^ T i[D) are finite constants. From ( 4.13Q we see that the second term 
appearing within square brackets in Ql.llj ) takes the following form for D — > 4: 



in 



2 

2 If 



(A-D) 



in* 



7 M In (an) + 7^ V + C M ({e M }, D) 



(4.16) 



Hence by combining ( 4.16 ) with ( |4.12| ) we find again that the terms proportional to 
(4 — D)~ cancel (as was the case in the previous two examples), and that we are left 
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with a /i-independent expression with the lattice spacing as the only regulator: 

-Q 2 

A M (p, p', m; a, {e M }, D = 4) = In (am) + 7m £> + Af^ (p, p', m) 

4 4 

+ ^tI^ItIq <Vrf? + t ^rlr V + (a) 

y,T.e=l r=l 

Here b^ T = b^ T (D = 4) is a finite constant, while a^ r e = a^T-e (D — > 4) = oo. Note 
that the non-covariant, {e^}- dependent terms, which destroy the renormalizability of 
the theory, vanish after averaging these expressions over all possible sets of {e^}. 

5 Conclusion 

In this paper we have studied in detail the fermion self energy, vacuum polarization 
tensor and vertex function for lattice QED in one-loop order, for the case where the 
derivative terms in the Dirac operator are replaced by one-sided lattice differences. 
Using the method of ref. we have systematically expanded the one-loop lattice 
Feynman integrals for small lattice spacing a up to O (a ). Although the small-a- 
expansion method makes use of the dimensional regularization of lattice Feynman 
integrals in intermediate steps, the in the limit D — » 4 divergent terms, possessing a 
covariant structure, were found to cancel. While the vacuum polarization tensor was 
found to have the correct continuum limit, the fermion self energy and vertex functions 
included non-covariant terms, which only vanished after averaging the expressions 
over all possible choices for the one-sided lattice differences. In the case of the vertex 
function these non-covariant terms were of O (a ) , while the fermion self energy also 
included a term of O (a -1 ). By discretizing the derivative in the Dirac operator using 
a fixed set of one-sided lattice differences, one is led to a non-renormalizable theory. 

Similar computations have been carried out by one of us (N.S.) for the case of 
Wilson fermions, where it is found that the contribution of O (a^ 1 ) to the fermion self 
energy has the structure of a mass term proportional to the Wilson parameter, which 
can be absorbed into a mass renormalization constant. No non-covariant contribu- 
tions were encountered in this case. 
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Appendix 



In all orders of the small-a-expansion we have studied, the infrared divergent contri- 
butions could be cast into the form 

+ " d D k 1 



M (D, n) 



(27T) 



D 



D 



4 £ sin 2 & 



(A.l) 



This expression is infrared divergent for n > 2 and D — > 4. The infrared divergent 



D 



part can be isolated as follows. Define q 2 = 4 J2 sin 2 -A Then 



P =i 



M(D,n) 



1) 



n-l 



1 



n-l oo 



U- (±\ [ dt e -i'r 
in-l)\J (2k) d W l( " 



o 



d D k 
7d 



dt i 



n-l 



-2t (l-cosfc^ 



(n-l)\l (2nf J 

Interchanging the t and ^-integration and setting s = 2t, one obtains 

1 x 



(A.2) 



M (D, n) 



2T in] 



ds s 



n-l 



l D 



e- s In s) 



(A.3) 



where In (s) is the modified Bessel function of zero order. 

Consider the case where n = 2. Then ( A.l ) is logarithmically divergent for 
D — > 4. The divergence arises from the large s behaviour of the integrand in (|A.3|) 
whose leading term, for n = 2 is given by (2tts) -D / 2 . We therefore decompose (|A.3|) 
for n = 2 as follows 



M(D,2) 



1 

4 



1 




D 00 ( 






J ds s 


e~ s In (s) 


+ \I dss { 


e- s In (s) 


D - 1 1 

(2tts)^ J 








4(2vr)-^ 



ds s 1 ? . 



(A.4) 
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The first two terms on the rhs are finite, and are denoted below by Af 2 . For D > 4 
this expression therefore takes the form 

Similar expressions can be obtained for n > 2. Such expressions are required for 
isolating the infrared divergent parts in higher orders of the small-a-expansion. For 
example the correction of O (a 2 ) to the vacuum polarization tensor calculated with 
fermionic action using one-sided lattice differences for the derivatives, requires the 
knowledge of M (D, 3), which can be shown to be given by 

where A/3 is a finite constant for D — > 4. 
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Figure 1: Diagrams contributing to the fermion self energy. 



P 





p — p 



a 



Figure 2: Diagrams contributing to the vacuum polarization tensor. 
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Figure 3: Diagrams contributing to the vertex function 
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